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1  Introduction. 


Let  {2’„,  n  >  1}  be  a  sequence  of  random  vectors  taking  values  in  such  that  Tnjn  — »  c 
in  probability  as  n  — >  oo.  When  the  closure  of  the  set  5  does  not  contain  c,  the  probability 
of  the  event  {Tn/n  E  5}  tends  to  0  and  oftentimes  at  an  exponential  rate.  Under  these 
circumstances  the  event  {Tn/n  E  5}  is  called  a  large  deviation  event.  Large  deviation 
theorems  provide  asymptotic  expressions  for  the  logarithm  of  probability  of  the  large 
deviation  events  while  strong  large  deviation  theorems  provide  asymptotic  expressions 
for  the  probability  of  the  large  deviation  events.  In  this  paper  we  will  be  mainly  concerned 
with  strong  large  deviation  results. 

In  the  univariate  case,  when  Tn  is  the  sum  of  n  independent  random  variables 
with  a  common  nonlattice  distribution  function  (d.f.)  F,  Cramer  (1938)  established  one 
of  the  earliest  strong  large  deviation  theorems  for  P(Tn/n  >  m).  The  lattice  case  was 
treated  by  Blackwell  and  Hodges  (1959).  Bahadur  and  Ranga  Rao  (1960)  generalized 
these  results  by  improving  the  remainder  terms,  both  for  nonlattice  and  lattice  d.f.  F 
under  the  assumption  of  finiteness  of  the  moment  generating  function  (m.g.f.)  of 
in  an  interval  around  the  origin.  When  the  m.g.f.  of  X\  is  finite  in  some  nondegener¬ 
ate  interval,  Petrov  (1965)  obtained  analogous  results  which  hold  uniformly  for  m  in  a 
compact  interval  with  zero  as  one  of  its  end  points.  Hoglund  (1979)  presented  a  unified 
formulation  of  the  results  of  Bahadur  and  Ranga  Rao  (1960)  and  Petrov  (1965),  both 
for  small  and  large  deviations.  He  gave  conditions  under  which  the  approximation  holds 
not  only  when  m  belongs  to  compact  sets  but  also  when  m  is  close  to  oo. 

A  k-dimensional  result  which  is  related  to  the  continuous  case  of  Hoglund  (1979) 
was  given  by  Borovkov  and  Rogozin  (1965).  These  authors  obtained  strong  large  de¬ 
viation  theorems  for  P(Tn/n  E  Sn)  under  minimum  conditions  on  the  measurable  sets 
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{5n,  n  >  1}.  Von  Bahr  (1967)  obtained  strong  large  deviation  theorems  for  measurable 
sets  {i>n,  n  >  1}  which  can  be  written  as  the  difference  of  two  convex  sets  considering 
two  cases.  In  the  first  case  Sn's  were  assumed  to  be  a  subset  of  a  sphere  with  its  center 
at  the  origin  and  in  the  other  case  Sn' s  were  assumed  to  be  contained  in  the  comple¬ 
ment  of  such  a  sphere.  In  a  recent  paper  Robinson,  et..  al.  (1990)  extended  the  results  of 
Hoglund  (1979).  They  considered  the  case  where  Tn  is  the  sum  of  n  independent  random 
vectors  in  Rk  with  the  first  ko  <  k  dimensions  being  lattice  with  span  1. 

In  this  paper  we  obtain  strong  large  deviation  theorems  for  arbitrary  sequences  of 
random  vectors,  under  simple  and  easily  verifiable  conditions  on  the  moment  generating 
functions.  These  results  for  arbitrary  random  vectors,  generalize  the  univariate  results 
of  Chaganty  and  Sethuraman  (1992)  and  those  of  Petrov  (1965).  The  key  to  the  strong 
large  deviation  result  is  a  local  limit  theorem  for  arbitrary  sequences  of  random  vectors, 
that  is,  a  theorem  on  the  convergence  of  pseudo-densities.  Before  obtaining  the  strong 
large  deviation  result,  we  will  prove  such  a  local  limit  theorem  under  mild  conditions  on 
the  characteristic  functions  (c.f.’s)  of  the  random  vectors. 

Let  {Tn,  n  >  1}  be  a  sequence  of  random  vectors  which  converges  in  distribution 
to  Y .  For  any  set  S ,  let  S°  be  the  interior  of  S  and  S  be  the  closure  of  S.  Let  5  be 
a  measurable  subset  of  Rk  such  that  0  <  y(S°)  =  y(S)  <  oo,  where  //  is  the  Lebesgue 
measure  on  Rk  and  bn  — >  oo  be  a  sequence  of  real  numbers.  Since  Yn  may  not  possess  a 
p.d.f.,  we  will  use  qn{y,  bn,  S)  =  bk  P(bn(Yn  —  y )  €  S)/f.i(S)  and  call  it  the  pseudo-density 
of  K.  The  convergence  of  qn{y,  bn,S )  to  the  probability  density  function  of  Y  will  be 
referred  to  as  a  local  limit  theorem.  In  this  paper  we  also  obtain  two  local  limit  theorems 
(Theorems  2.1  and  2.2)  for  random  vectors  Yn  under  different  conditions  on  the  c.f.  of 
Yn.  A  key  result  concerning  asymptotic  expressions  for  the  Laplace  transform  on  the  set 
Yn  >  0  is  proved  in  Theorem  2.4.  This  result  is  used  later  in  Section  3. 
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As  mentioned  earlier,  these  local  limit  theorems  are  used  to  obtain  a  strong  devi¬ 
ation  result  (Theorem  3.4)  for  arbitrary  random  vectors.  As  an  application  of  this  large 
deviation  result,  we  obtain  expressions  for  the  tail  probabilities  of  the  k  dimensional 
multivariate  F-distribution  in  Section  4. 


2  Multidimensional  Local  Limit  Theorems. 


Let  {V„,  n  >  1}  be  a  sequence  of  random  vectors  taking  values  in  Rk ,  which  converges 
to  a  random  vector  Y  in  distribution.  Let  5  be  a  measurable  subset  of  Rk  satisfying 
0  <  /i(5°)  =  fi(S)  <  oo,  where  y.  is  the  Lebesgue  measure  on  Rk.  Let  bn  — »  oo  be  a 
sequence  of  real  numbers.  Define 

q4y,b„,S)  =  JLp[bn(Yn-y)€S)  (2-1) 

for  y  in  Rk.  Notice  that  qn{y,  K,  S)  is  the  p.d.f.  of  +  Zn  where  Zn  is  distributed 
uniformly  on  —S/bn.  We  will  refer  to  it  as  the  pseudo-density  of  Yn  and  study  its  limiting 
properties,  since  Vn  itself  may  not  possess  a  p.d.f.  Let  {y„,  n  >  1}  be  a  sequence  of 
random  vectors  in  Rk  such  that  yn  — >  y‘.  We  will  refer  the  convergence  of  yn(yn;  bn,  S ) 
to  the  p.d.f.  of  Y  at  y’  as  a  local  limit  theorem.  We  will  use  the  following  notations  to 
denote  norms,  vector  products,  etc.  Let  t  =  (fj,. . .  ,f*),  s  =  (sl5...,s*)  be  vectors  in 
Rk .  Define 


<t,S>  —  tiSi  +  •  •  •  -f  tkSk 

11*11  =  max  |fj! 

1*1  =  |*i|  H - +  \tk\ 

ts  = 

t  >  s  if  ti  >  Si  V  i.  (2-2) 
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We  first  obtain  a  local  limit  theorem  under  a  strong  condition.  This  is  given  in 
Theorem  2.1  below.  This  strong  condition  is  relaxed  and  a  more  useful  result  is  given  in 
Theorem  2.2. 

THEOREM  2.1  Let  {y'n,7i  >  1}  be  a  sequence  of  random  vectors  which  converges  weakly 
to  a  random  vector  Y.  Let  fn(t)  and  f(t )  be  the  c.f.  's  ofYn  and  Y  respectively.  Suppose 
that,  there  exists  an  integrable  function  f’(t)  sxich  that 

sup|/„(OI  <  no  (2-3) 

n 

for  all  t.  Then  Yn  possesses  a  bounded  and  continuous  p.d.f.  fn{y),  V  also  possesses  a 
bounded  and  continuous  p.d.f.  f{y),  and  fn(yn)  converges  f(ym)  if  yn  — ►  y' ■ 

Proof.  Condition  (2-3)  implies  that  the  c.f.’s  fn(t)  and  f(t)  are  integrToie.  Hence 
both  Yn  and  Y  possess  a  bounded  and  continuous  p.d.f.'s.  The  inversion  formula  and 
the  dominated  convergence  theorem  show  that  fn{yn)  converges  to  f(ym)  if  yn  ->  y‘  as 
n  — >  oo.  0 

Condition  (2-3)  is  too  strong  to  be  useful  in  most  situations.  We  show  in  The¬ 
orem  2.2  below  that  appropriate  bounds  on  the  c.f.  fn{t)  on  increasing  sequences  of 
bounded  intervals  are  sufficient  to  obtain  results  similar  to  Theorem  2.1. 

THEOREM  2.2  Let  {Yn,  n  >  1}  be  a  sequence  of  random  vectors  which  converges  to  a 
random  vector  Y  in  distribution.  Let  fn(t)  be  the  c.f.  of  Yn  for  n  >  1  and  let  f(t)  be 
the  c.f.  of  Y .  Let  {/?„}  and  {&„}  be  sequences  of  real  numbers  such  that  j3n  — >  oo  and 
bn  — >  oo.  Suppose  that  there  exists  an  integrable  function  fm(t)  such  that  for  t  €  Rk . 

sup  |/»(«)|  /(||*||  <  A.)  <  /*(*)  (2-4) 
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and 


MA)  ^  SUp  |/„(<)|  =  0(71 

0n<||l||<A!>„  \bn/ 


(2-5) 


for  each  A  >  0,  where  the  above  supremum  is  defined  to  be  zero  if  {t  :  fin  <  ||l||  <  \bn }  is 
empty.  Then  the  random  vector  Y  possesses  a  bounded  and  continuous  p.d.f.  f{y).  Let 
qn(y,  bn,  S)  be  the  pseudo-density  function  of  Yn  as  defined  in  (2-1).  Then  there  exists 
a  finite  constant  M  and  an  integer  ns  depending  on  S.  such  that 


sup  [qn(y\  6n,  5)]  <  M 


(2-6) 


for  n  >ns.  Furthermore,  if  yn  — >  y *  then 


qn(yn ;  S)  f{ym) 


(2-7) 


as  n  — >  00. 

Proof:  Since  fn{t)  — >  /(<)  pointwise  and  fin  — >  00,  condition  (2-4)  implies  that  f(t)  is 
bounded  by  /*(<)■  Hence  V  possesses  a  bounded  and  continuous  p.d.f.  f(y).  Suppose 
/ 3n/bn  is  bounded.  Since  b^6n( A)  — »  0,  for  each  A  >  0,  we  can  find  a  sequence  {A„} 
satisfying 


An  — >  00  and  A^6^6»T1( A„)  ->  0  (2-8) 

as  n  — »  00.  If  fin/bn  —y  00,  put  An  =  ^n/6n  and  in  this  case  also  (2-8)  is  satisfied  because 

def 

0n(An)  =  0  for  large  n.  For  simplicity  of  notation  set  6n  —  6n( An).  Let  Un  be  the 
uniform  distribution  on  the  set  —S/bn  and  un,  u„  be  the  p.d.f.  and  c.f.  corresponding 
to  Un.  We  also  introduce  another  distribution  function  (d.f.)  with  p.d.f.  un  and  c.f. 
vn,  to  obtain  the  important  identity  (2-11): 
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v  (x)  -  TT  [sin(An6njj/2)l  <  X  <  oo  j  =  1  A-  (9.g) 

tn(l)  "  (2x)* /j  [  (Wi/2)  J  ’  J 


»»n(<)  = 


n  i- 


if  IUII  <  A nbn 


otherwise. 


(2-10) 


Let  Fn  be  the  d.f.  of  ln,  and  let  Qn  =  Fn  *  Un,  Mn  =  Qu  *  Vn  where  *  denotes  the 
convolution  operation.  Notice  that  qn(y,  bn ,  S)  defined  in  (2-1)  is  the  p.d.f.  of  Qn.  Let 
mn(y)  be  the  p.d.f.  of  Mn.  The  c.f.  ra„(i)  of  Mn,  which  is  equal  to  fn{t)un{t)vn(L ), 
vanishes  outside  the  rectangle  {t  e  Rk  :  \\t\\  <  An6„).  The  inversion  theorem  yields  the 
following  identity: 


exp(-i  <t,  y  >)  mn(t)  dti  . . .  dtk 


ln  (y)  =  J7TTZ  [  "■  f  exp  (-i  <  t,  y  >)  mn(t)  d 

(27T )*  J  j\\t\\<\nb„ 

/OO  TOO 

•••/  qn(y-x;bn,S)vn(x)dxl...dxk 

-oo  J — OO 


=  [  •••  /  P(bn{Yn-y  +  x)e  S)vn(x)dx1...dxk.  (2-11) 

)  J -co  —  oo 

Relation  (2-11)  is  the  starting  point  of  the  main  part  of  this  proof  and  it  relates 
qn(y;  bn,  S )  to  the  integrable  c.f.  m„(<).  Let  yn  — ►  y*.  We  first  show  that  mn(yn) 
converges  to  /(y*)  and  then  obtain  lower  and  upper  bounds  for  mn(y„ )  which  depend 
on  qn{yn',  bn,  S).  This  will  then  establish  (2-6)  and  (2-7).  Using  (2-8)  and  (2-11)  we  get 


in(y)  —  /ri  \k  f  '  '  '  [  exP(  *  ^  F  y  )  ^n(A)  dt\  . . .  dtk 
(27 r)*  J  -/||t||</?n 


(Stt)*  /  Ipn 


<|h||<A„tn 


exp(— i  <  t,  y  >)  ihn(t)  dt\ . 


1  r03 


/co  r  oo  A  *  hr  B 
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(2-12) 


9  fee*  r  oo 

<  ■■■/  nt)dh...jtk. 

(2n)K  J- oo  J- co 

for  ?i  >  n0,  where  no  is  independent  of  y.  From  condition  (2-4),  (2-8),  (2-12)  and  the 
dominated  convergence  theorem  and  the  inversion  formula  we  get 

1  r oo  roo 

mn(yn)  — *  -  /  •••/  exp  (-i.  <  t,  y'  >)  f  {t)  diy  . . .  dtk  =  f(ym).  (2-13) 

(Z7T  j  J —co  J —oo 

Let  rj  >  0.  Let  S'(.t,  rj)  —  {y  :  j|y  —  xjj  <  ?/}  be  a  ball  of  radius  ?/  (with  norm  ||  •  j|) 
centered  at  x.  Let  Sv  =  {x  :  S(x,  tj)  C  S}  and  S’1  —  {y  :  ||x  —  y||  <  y  for  some  x  €  5). 
Since  we  have  assumed  that  y(S°)  =  y(S)  we  can  find  y(=  ys)  >  0  such  that 

fi{Sr,)  >  0  and  [/x(5T?)/y(S)]  <  2.  (2-14) 

Note  that  y  €  Sv  implies  y  +  x  €  S(y,  y)  C  S  if  ||x||  <  y.  From  this,  we  get  a  lower 
bound  for  mn(y)  as  follows: 


mn{y)  >  /•••/  P(bn(Yn-y +  x)  e  S)vn{x)dx1...dxk. 

Ill  )  J  A\x\\<r)/bn 


f\\x\\<7)/bn 


KS) 


> 


KS) 

K 

y(S) 


P(bn(Yn  -  yn)  e  5,) 


1  - 


|x||<Tj/tn 
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«»(*)  dxj  .  ..dxk 


T7\nTj 


(2-15) 


Using  (2-12),  (2-16)  and  (2-8)  we  get 


V(S)  [  7rAn7/J 

1  roo  f  oo  \kh^0 

<  WL-L™dt'-dt>+W 

2  /"CO  fOO 

<  75^5  /  •••/  (2-16) 

j  J —co  J — co 


/ 


for  sufficiently  large  n.  By  replacing  S  by  S'1  and  using  (2-15)  and  the  fact  that  S  C 
we  get 


< 


J?Lp(bn(Yn-yn)£Sv) 

niS")  2 


7rA„7/ 


')  rc o  r  co 

3* •  ••  y_  r[t)dtl...dtk 


fl(S)  (2*y 

4  roo  roo 

5  (2 7? L-L™ 

Since  An  — >  co  as  ??  — >  oo  we  can  find  an  integer  ns  so  that 


1  bt 


sup 

y 


KS) 


3  ^(fcn(>n  -  y»)  6  S„) 


<  M 


(2-17; 


(2-18) 


for  n  >  ns,  where 


M  =  TTYk  r  •••  [°°  nt)dU---dtk-  (2-19) 

(J7T)  J —oo  J—  oo 

This  proves  assertion  (2-6).  Note  that  y  €  S  implies  that  y  —  x  €  Sv  for  ||x||  <  rj. 
Therefore  for  n  >  ns  an  upper  bound  for  7nn(y)  is  given  by 


mn{y)  = 

< 


< 


bt 


KS) 


A  roo  roo 

3  /  ■■■  P(bn(Yn  -  y  +  x)  €  S)  vn{x)  dxx  . .  ,dxk 

O  )  J  —  OO  J  —  OO 

un(x)  dx!  . . .  dxk 


y{S) 

| *•)/-/ 


+71/ 


J-k 


l|r||<J}/6n 
vn(x)  dxi  . . .  dxk 


-—j  P(bn(Yn  -  yn)  €  5”)  +  M 


1  - 


1  - 


7rA„77 


(2-20) 


Thus,  from  (2-14),  (2-16)  and  (2-21)  we  get  that 
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limsup  P(bn{Yn  -  ;/„)  £  Sv) 

n  lfi(b) 

<  /(y*)  <  lim  inf  P{b,t(Yn  - 

L  n  Z'P) 

Replacing  5  by  Sn  in  the  l.h.s.  and  S  by  Sv  in  the  r.h.s.  and  using  the  relations  S'  C  (5’’% 
and  (Srj)71  C  S  we  get  that 

limsup  -73-7  P(bn{Yn  -  y„)  £  S) 

n  [p(STI) 

<  /(y')  <  lim  inf  P(bn(yn  _  yn)  €  <?)j  .  (2-22) 

71  IviSr,) 

Letting  77  — >  0  and  using  the  fact  fi(S°)  =  p(S)  we  get  the  assertion  (2-7).  <C> 

We  have  replaced  the  strong  condition  (2-3)  of  Theorem  2.1  with  a  weaker  but  more 
comlicated  looking  condition  (2-4)  in  Theorem  2.2.  The  next  theorem  gives  a  convenient 
condition  to  verify  condition  (2-4). 

THEOREM  2.3  Let  {V’n,  n  >  1}  be  a  sequence  of  random  vectors  taking  values  in  Rk , 
with  c.f.’s  { fn (0 }  •  Let  {dn}  be  a  sequence  of  real  numbers  such  that  dn  — >  00.  Assume 
that  there  exists  6  >  0  such  that  gn{t)  =  dff2  log  |/n(c/nt)|  is  differentiable  in  an  open 
neighborhood  U  =  {t  :  ||t||  <  6}  of  the  origin  for  all  n  >  1.  Suppose  that  there  exists 
a  >  0  such  that  for  t  £  U , 

t'V2gn{t)t  <  —at't  (2-23) 

for  all  n  >  1.  Then  condition  (2-4)  of  Theorem  2.2  is  satisfied  with  8n  =  6dn. 
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Proof:  An  application  of  Taylor's  theorem  yields  for  t  £  U. 

n\  io\  i  ^  j  V7  /n\  ^  i  ^  ^  9n{>  n  )t 

9n{t)  =  <  t,Vgn(0)  >  + - - - 

/'V2gn(rn)<  <  ai't 
~  2  _  2 

where  rn  £  £/.  Therefore  for  t/dn  £  £L  we  have, 

<  .  ^  oT/ 

9n(t/dn)  <  -JJT- 

^”71 

Let  /?„  =  <Mn.  Thus 

LMOIWII  <  M  <  exp(4(S„((/4))) 

<  exp(—ai't/2) 

which  is  an  integrable  function.  This  completes  the  proof  of  the  theorem.  <0 

The  next  theorem  obtains  the  limit  of  a  function  related  to  the  Laplace  transform 
of  the  positive  part  of  Yn  when  (2-6)  and  (2-7)  hold.  It  plays  an  important  role  in  the 
proofs  of  the  strong  large  deviation  theorem  of  Section  3. 

THEOREM  2.4  Let  { Tn ,  u  >  1}  be  a  sequence  of  random  vectors  converging  weakly  to 
Y .  Let  {6n}  be  a  sequence  of  real  numbers  such  that  bn  — ►  oo.  Let  S  be  a  measurable 
subset  of  Rk  such  that  0  <  ft(S°)  =  //(5)  <  oo.  Let  qn(y: ;  bn .  S)  be  as  defined  in  (2-1). 
Assume  that  qn(y,  bn ,  S )  satisfies  (2-6)  and  (2-7).  Then 

6^[exp(-6nX:Vn,)/(V'n>0)]  -  /( 0)  (2-27) 


Proof:  Let  h  >  0  be  given.  Consider 

k 

ln  =  £[exp(-6n^yni)  /(}„  >  0)] 


(2-24) 

(2-25) 

(2-26) 
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=  EE -E  e 

<i=i<s=i  4=1 

Let  ynit  =  (2 /,  —  l)i/2in.  Then 

h 


oxp(-4„f  '«,)  n  1 1 — LIE  <  >v  <  ,i*' 

i=i  i=i 


in 


in 


(2-28) 


CO  CO  CO 


EE -E  £ 

<1=1  L=i  <(.=1 


•  k  ( “i  /i 

e^P(  in  ^  ^«*)  IT  ^  I  7  —  ^  ni  JJnl,  8  T 


i=l 


i  =  l 


(2-29) 


,2h„  2in< 

Let  ?n  =  m/,  =  [1/i2]  and  5'/,  =  [-i/2,  i/2]\  /  =  (/i,. . .  ,/*)  and  j/„/  =  (j/n<,, - ?yn/fc ) - 

We  now  get  lower  and  upper  bounds  for  In  as  follows: 

mm  m  k 


i.  >  EE  '■■£e*p(-'>£'.)r(-;^<n,.-9.,<.  <2-.v; 

<1=1  <2  =  1  4  =  1  1  =  1  ^ 


I  k  m  m  m  /  k  \ 

jfEE'-Z  exp  (  —  h  EM  <7n(?/n<i  ini  ‘8<i  ) 


(2-30) 


n  <1=1  <2=1  4  =  1 


1=1 


and 


/A-  oo  oo  oo  /  fc  \ 

4  <  TT  H  I]  •  •  ■  5Z  eXP  -/?  -  !)  9n(yn<;  in,  5fc) 

"n  l^lijrl  4=1  V  t=l  / 

1  fc  m  rn  m  /  k  \ 

=  TT  S  £  •  •  •  5T  eXP  _/t  -  !)  )  (?/„<;  in.  50 

«  /]  =1  <2  =  1  4  =  1  \  1  =  1  / 

+^^^<hJF  exp  (-/z  Jy(/’  ~ 1})  7T,(j/n';  6n’  5h)' 

Since  r/„;  — >  0  for  each  /  we  have  qn[ynl\  in,  5’/i)  — >  /(0)  as  n  — *  oo.  Thus,  we  get 

i  fc 


(2-31) 


lira  inf  (i*/n)  >  /(0)/<; 


=  /(0) 


£exp(— Ih) 

Ll=l 

i(exp(  —  h)  —  exp(  —  h(m  +  1))) 


1  k 


1  —  exp(  — i) 


and 


limsup(*;/„)  <  /(0)[/<XX,exp(-(l-  1  )/•)]*  +  iir^S»=rS£„+1exp(-(/-  1)4) 
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=  AO!  [M1,rgi:Ti  “1 '  +  *'*  (rsfed  ™i>< ->»*)■  (2:«) 

From  (2-33)  and  (2-34).  letting  h  — >  0  and  noting  that  ?n  =  []///2].  we  get 

lim(^;/n)  =  /(0).  (2-34) 

This  completes  the  proof  of  the  theorem.  <C> 

We  will  need  the  following  extension  of  Theorem  2.4  in  the  next  section. 


THEOREM  2.5  Let  {V^,  r?.  >  1}  be  a  sequence  of  random  vectors  in  Iik .  Let  fn(t)  be 
the  c.f.  of  Yn  for  n  >  1.  Let  {Sn}  and  { bn }  be  sequences  of  real  numbers  such  that 
(3n  — ►  oo  and  bn  —*  oo.  Assume  that  {>'),,  n  >  1)  satisfies  conditions  (2-4)  and  (2-5) 
of  Theorem  2.2.  Let  {Pn,  n  >  1}  be  a  sequence  of  real  positive  definite  matrices  whose 
eigenvalues  are  bounded  above  and  bounded  away  from  0.  Let  P)J2  be  the  unique  positive 
definite  square-root  of  Pn.  Suppose  that  P^2Yn  converges  in  distribution  to  a  random 
vector  V'.  Then 

IPn^'^bl  E[exp(-bnJ2Yni)  /(};  >  0)]  ->  f(0)  (2-35) 

t=i 

as  n  — >  oo,  where  f  is  the  p.d.f.  of  Y . 


Proof:  We  will  show  that  every  subsequence  has  a  further  subsequence  such  that  (2- 
36)  holds  for  the  subsequence.  Consider  a  subsequence  and  denote  it,  for  simplicity,  as 
the  sequence  {?r}  itself.  Since  the  elements  of  P„  are  uniformly  bounded,  we  can  find 
a  further  subsequence,  which  again  for  convenience  is  labelled  as  {77},  and  a  positive 
definite  matrix  P  such  that  Pn  —y  P  as  n  — *•  00.  Let  Q  =  P~XP ,  then  Yn  converges  to 
Q  Y  in  distribution.  Note  that  the  p.d.f.  of  QY  at  zero  equals  |(5|/(0),  where  /  is  the 
p.d.f.  of  Y.  The  conclusion  (2-36)  now  follows  from  Theorems  2.2  and  2.4.  0 
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3  Strong  Large  Deviation  Theorems. 


Let  {Tn,  v  >  1}  he  a  sequence  of  random  vectors  in  Rk .  Assume  that  the  m.g.f.  <?„(-)  = 
£’[exp(<  z,  Tn  >)]  is  holomorphic  and  nonvanishing  in  flfc(c),  where  fi(c)  =  {~  E  C  : 
|c|  <  c }  for  some  c  >  0  and  C  is  the  set  of  complex  numbers.  Let  {an}  be  a  sequence  of 
real  numbers.  Then 

n>„{z)  =  —  log0„(*)  (3-1) 

on 

is  a  well  defined  holomorpnic  function  on  Clk(c).  Let  VV’n(*)  =  (Diij.'n(z)'  . . . ,  DkU'n{~)) 
be  the  vector  of  first  order  partial  derivatives  and  V2il'n(z)  denote  the  matrix  of  second 
order  partial  derivatives,  that  is, 

V2fc(~-)  =  (D.jMz)).  (3-2) 

The  determinant  of  the  matrix  V2i/>n(~)  is  denoted  by  |V2i/>n(~)l-  For  u  in  Rk .  let 

7n(u)  =  SUp[<  S,  U  >  -li'n(s)}.  (3-3) 

se/fc 

Let  {mn}  be  a  bounded  sequence  of  vectors  in  Rk  such  that  there  exists  a  sequence  of 
vectors  {t„}  satisfying 

Vt/>n(Tn)  =  mn  and  0<d<rnt<b<c  for  all  1  <  i  <  k.  n  >  1.  (3-4) 

Under  these  conditions  we  can  see  that  7„(m„)  =<  m„,  rn  >  —d>n(Tn).  Let  Kn  be 
the  d.f.  of  Tn.  We  will  use  the  left  continuous  version  of  the  distribution  function  which 
will  enable  us  to  write  the  identities  in  (3-6).  Let 

dlUv)  =  eX'P>^  f’  dl\„{y),  yeR*  (3-5) 

r*n  (  Tn  ) 
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and  T’  be  random  variable  with  d.f.  Hn{y).  Let  T'n  =  (7'*  —  anmn)/ y/a^  and  U„  = 

•  •  • ,  Ynk)  "here  V’m  =  rni7,t'1.  Using  these  new  random  variables  we  can  write 

P{Tn  >  anmn)  =  f  dl\n(y) 

J{y£Rk:y>anmn} 

/  <  V  •)  i  n  >  ~\~&n  V-yn  (  1  n  )]  dlln{y) 

J  {y€Rk-y>anmn} 

=  exp(a„^B(rn))  7?[exp(-v^'  <  rn,  T’  >)  I(I\ '  >  anm„)] 

=  exp(-a„7n(mn))  £[exp(-  ^  rni  T'ni)  I (T'n  >0)] 

t  =  l 

k 

=  exp(-an7n(mn))  £[exp(-6n  ]T  *m)  I{Yn  >  0)]  (3-6) 

i=i 

where  bn  =  a ~  and  Yn  is  as  defined  before.  The  next  lemma  obtains  some  important 
properties  of  the  sequence  { Yn ,  n  >  1}. 

LEMMA  3.1  Let  {Tn,77  >  1}  be  a  sequence  of  random  vectors  in  Rk .  Let  {m„ }  be 
a  bounded  sequence  of  vectors  in  Rk  such  that  there  exist  a  sequence  of  vectors  {rn} 
satisfying  (3-4)-  Let  the  random  vectors  {Yn,n  >  1}  6e  as  defined  above.  Assume  that 
{7"n,  n  >1}  satisfies  the  following  conditions: 

(A)  There  exists  0  <  /?  <  oo  such  that  |V’n(-)|  <  3  for  all  n  >  1,  z  £  fU(c). 

(B)  There  exists  a  >  0  such  that  the  eigenvalues  of  V2^’n(rn)  are  bounded  below  by  a 
for  all  n  >  1 . 


Let  Dn  be  the  diagonal  matrix  with  rni  as  the  ith  diagonal  element.  Then,  Y'n  is  asymptot¬ 
ically  multivariate  normal  with  mean  0  and  variance-covariance  matrix  Dn  V2t/>„(rn)  Dn. 


Let  fn(t)  be  the  c.f.  oJYn.  Then  there  exists  6  >  0  and  a ,  >  0,  independent  of  n,  such 
that 


sup  \fn(t)\ I{\\t\\  <  hy/af)  <  cx p(  ex\t' t).  (3-7) 


Proof:  The  c.f.  of  T'n  is  given  by 


<0(0  =  exp (-iy/eu  <  i ,  mn  >) 


tO ( rn  T  it  j yjet n ) 

Qn  ( Tu ) 


(3-8) 


Since  tpn{~)  is  a  holomorphic  function  in  f0(c)>  for  )f|  <  (c—  b)j 2  we  have 


log  </n(0  =  <  0  mn  >  +«n[^n(T„  +  ? 


<'V2V-n(r„)t  ,  _  ,  ,  *  N 

0  "t"  ■*?n(^Vi  ”f“  2  t - ) 

\I CL  n 


tf’n  ( T’n  )  ] 


(3-9) 


where  /?n(rn  -f  =  YffyL?,Yi{a  ■.  |a|=j}  aLnH*0o>  By  Cauchy’s  theorem  and  condition  (A) 
(see  (3.7)  of  Chaganty  and  Sethuraman  (1986))  we  get  the  following  bounds 


|Aj0n('Tn)|  <  t)2,  |flin)|  <  ^c_^|0|-  (3*10) 

Thus  the  elements  of  the  matrix  V2t/>„(Tn)  are  bounded  uniformly  in  n.  Now  for  |<|  < 
(c  —  6)/2  and  for  all  n  >  1,  we  have 


Hence 


|tfn(Tn  +  *01  < 

< 


CO 


k 


E  E  ni'.r 

J=3  [o  ,  I i=i 


<*n  j  -^n  (  *"n  T  * 


oo 

< 


j=3 


liP  jk 

(c  -  02  <0/2 


(3-11) 
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< 


(3-12) 


as  ??  — *  oo.  From  (3-9)  and  (3-12)  we  have  that  T'u  is  asymptotically  multivariate  normal 
with  mean  0  and  covariance  matrix  V2t/>„(t„).  Let  Dn  be  the  diagonal  matrix  with  rnt 
as  the  ?'th  diagonal  element.  Note  that  Yn  =  Dn T'n  and  V’n  is  asymptotically  multivariate 
normal  with  mean  0  and  covariance  matrix  DnV2ii'n(Tu)  D„.  We  now  proceed  to  show 
(3-7).  Let 


9n{t)  =  (1^  l0S  l$»(v^OI 

=  Real  {u>n(Tn  +  it)  —  Vn{Tn))-  (3-13) 


We  find  that. 


t'V2gn(t)t  =  —  i'Real  (V2?/!n(Tn  +  it))  t 

=  —t1  V2i/’n(rn)  t  +  Real(7?T,(rtl  +  it)).  (3-14) 

From  (3-11)  we  have  |i2n(rn  +  i<)|/|/|2  converges  to  0  as  |<|  — >  0  uniformly  in  n.  Therefore 
we  can  find  a  61  >  0  such  that  for  ||£||  <  6i, 

\Rn{rn  +  it)\<  -^\t\2  (3-15) 

Frorr  (3-14),  (3-15)  and  condition  (B)  we  have 

t'V2gn(t)t  <  —at't  +  | t't  =  t't  (3-16) 

for  ll^ll  <  8\.  This  shows  that  the  c.f.  gn{i)  satisfies  condition  (2-24)  of  Theorem  2.3. 
Hence 


sup  |£„(0I  /(INI  <  hy/chl)  <  exp(— qT  t/2).  (3-17) 
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The  c.f.  fn{t)  of  Yn  is  given  by  fn(t)  =  gn(Dnt).  From  (3-1)  we  get  that  d\\t\\  <  HD,,  /||  < 
a0|NI  and  t Pi' t  <  t'DnDnt  <  a^l'l.  Let  6  =  b\/«o  and  =  ad2/'!.  Then  from  (3-17) 
we  get 

sup  |  fn (01  /(INI  <  ^  exp(-aN0-  (3-18) 

n 

This  completes  the  proof  of  the  Lemma.  0 


LEMMA  3.2  Let  {T„,  n  >  1 }  be  a  sequence  of  random  vectors  in  Rk .  Let  |r/?„}  be 
a  bounded  sequence  of  vectors  in  Rk  such  that  there  exists  a  sequence  of  vectors  {r„} 
satisfying  (3-f)-  Let  an  —>  oo  be  a  sequence  of  real  numbers  and  let  bn  =  y/af.  Suppose 
that  Tn  satisfies  conditions  (A)  and  (B)  and  the  following  condition: 


(C)  There  exists  S0  >  0  such  that 


sup 


<?n(rn  +  t'O 
<Pn{Tn) 


(3-19) 


for  any  given  6  and  X  such  that  0  <  6  <  S0  <  A.  Let  Yn  be  as  defined  in  (3-6).  Then 

k  k 


n  r"*  lVVn(T„)|1/2  bkn  £[eXp(-6„  In.)  /(Fn  >  0)] 


1 


i=l 


i  =  l 


(27t)*/2 


(3-20) 


Proof:  We  will  prove  this  lemma  by  verifying  the  conditions  of  Theorem  2.5  for  the 
sequence  {Fn}-  From  Lemma  3.1,  Yn  is  asymptotically  multivariate  normal  with  mean 
0  and  covariance  matrix  Dn  V20n(rn)  Dn.  From  the  same  lemma,  there  exists  a  6  >  0 
such  that  (2-24)  holds  with  /?„  =  6bn.  Using  condition  (C)  we  get  that  for  fixed  A  >  0, 

sup  I  .fn  (01  =  SUP  \9n{Dnt)\ 

6i„<||t||<Afc„  Sb„<Htll<Ab„ 

<  sup  \gn(t)\ 

6dbn<\\t\\<\bbn 


17 


(3-21) 


SU|> 


On  (  *7t  pit) 


On{Tn) 


This  verifies  that  t lie  c.f.  of  Vn  also  satisfies  condition  (2-25)  of  Theorem  2.2.  From  (3-4). 
(3-10)  and  condition  (B)  we  can  see  that  the  eigenvalues  of  Pn  are  bounded  above  and 
away  from  zero  below.  Thus  Pn  satisfies  the  conditions  in  Theorem  2.5.  The  assertion 
(3-20)  now  follows  from  (2-36)  of  Theorem  2.5.  0 


REMARK  3.3  The  conclusion  (3-20)  of  Lemma  3.2  can  also  be  rewritten  in  the  form 

k  1 


bkn  E[exp(  — 6n  £>;,)/(}•'„  >0)] 


(3-22 


(2tt )fc/2  nf=1  rni  |V2t'n(r„)|1/2 
where  the  symbol  ~  means  that  the  ratio  of  the  two  terms  above  tends  to  1  as  n  — ►  oc. 


We  are  now  in  a  position  to  state  the  main  theorem  of  this  section. 

THEOREM  3.4  Let  { Tn ,  n  >  1}  be  a  sequence  of  random,  vectors  in  Rk .  Let  {mn}  be 
a  bounded  sequence  of  vectors  in  Rk  such  that  there  exists  a  sequence  of  vectors  {rn| 
satisfying  (3-f).  Let  an  — >  oo  be  a  sequence  of  real  numbers.  Suppose  that  {Tn,  n  >  1} 
satisfies  conditions  (A),  (B)  and  (C).  Then 


'T 

P  (  —  >  m, 


exp(— a„7„(mn)) 


(2tt an)kl*  n,=i  rni  |VVn(rn)|1/2 


(3-23) 


Proof:  Let  Yn  be  a  function  of  Tn  as  defined  in  (3-6).  The  identity  (3-6)  states  that 

k 

P{Tn  >  anmn)  =  exp(-an7n(mn))  E[exp(-bn  ^  Fni)  I(Yn  >  0)]  (3-24) 

t=i 

Conditions  (A),  (B)  and  (C)  imply  that  {}j,}  satisfies  the  conditions  of  Lemma  3.2. 
The  conclusion  (3-23)  follows  from  (3-24)  and  (3-22).  0 
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4  Applications. 


In  this  section  we  will  demonstrate  the  applicability  of  the  theorems  in  Section  3  with 
an  example.  This  example  is  concerned  with  approximating  the  probability  of  the  tail 
areas  of  the  multivariate  F-distribution. 

Let  G\, C'kn  and  H„  be  independent  and  distributed  as  chi-square  with  ?/  A)  At 

and  n  degrees  of  freedom  respectively,  where  A),..., A*  are  fixed  positive  integers.  Let 
Gn  =  (Gi„, . . . ,  G\n).  The  distribution  of  the  random  vector  An  =  G„/lTn  is  known  as 
the  multivariate  F-distribution  with  (uAi,. . . . ?;A k,  n)  degrees  of  freedom.  Fix  a  vector 
r0  =  (roi, . . .  ,r0jt).  We  will  use  Theorem  3.4  to  estimate  the  tail  probability  P(Xn  >  r0), 
which  can  be  rewritten  as 


P(Xn>r0 )  =  P(Gn-r0Wn>  0) 
=  P(Tn>  0) 

where  Tn  =  G„  —  r0Wn.  The  m.g.f.  of  Tn  is  given  by 


(4-1) 


<t>n(Z) 


-n(nW) 


A,n/2 


n/2 


,l  +  2Ef=i^0L 


C  e  n* 


(4-2) 


where  H  =  {z  €  C  :  Real(z)  <  1/2}.  We  will  verify  that  {Tn}  satisfies  the  conditions 
(A),  (B)  and  (C)  with  an  =  n  and  mn  =  0  for  n  >  1.  Note  that  t/.'n(r)  =  A  log  ^n(r)  is 
independent  of  rr  and  is  given  by 


and  *l’(z)  is  bounded  for  Qk{c)  for  any  0  <  c  <  1/2.  Thus  {Tn}  satisfies  condition  (A). 
Since  mn  =  0  and  tpn  =  %l>  are  independent  of  n,  we  have  rnT  satisfying  (3-4)  do  not 
depend  on  n  and  are  equal  to  r0,  where  the  ?th  component  of  t0  is  given  by 


To,  = 


1 
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A,(l  +  kf) 
r0,(l  +  k\) 


1  <  *  <  k. 


Assume  that  the  vector  r0  —  (r01, . . . ,  r0jt)  satisfies 


(1-4) 


To,  >  A, 

1  T  k  v  1  +  A;  A 


for  all  1  <  i  <  k 


(4-5) 


where  f  =  J2ro,/k  and  A  =  A,//:.  From  (4-4)  and  (4-5)  we  can  verify  that  there  exists 

d  >  0  and  0  <  b  <  1/2  such  that  d  <  r0,  <  6  for  all  1  <  i  <  k.  Thus  r„  =  r0  satisfies 
(3-4).  With  some  simple  calculations  we  can  show  that  the  determinant  of  V2V'("o)  is 
given  by 


|vV(To)l = n  (if )  (iTi?f  (1+t'X)  <«> 

showing  that  the  eigenvalues  of  V2^’(r0)  are  bounded  below  away  from  zero.  Thus  { Tn } 
satisfies  condition  (B).  Using  (4-2)  we  can  see  that  for  any  given  6  >  0  and  A  >6, 


sup 

{t:6<||(||<A} 


<t>n(To  +  it) 

4>n{T0) 


<  sup 


Mto  ±  if) 

Mto) 


(4-7) 


where  <f)0 (z)  =  nf=i(l  —  2 z,)-*’/2  is  the  joint  m.g.f.  of  k  independent  chi-square  random 
variables.  Since  ^>o(To  +  t i ) / <^>o ( "^o )  is  the  c.f.  of  a  nonlattice  random  vector,  there  exists 
a  T]  >  0  such  that 
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sup 

{«  :fi<INI<M 


(4-8) 


Mto  +  *0 


4>o{t0) 


<  e 


-v 


From  (4-7)  and  (4-8)  we  get  that 


sup 

{t:6<||t||<A} 


<?o(to  +  f-0  _  f  1  \ 
<Po(t o)  °\nkl2) 


(4-9) 


Thus  we  have  verified  that  {Tn}  satisfies  condition  (C).  We  now  can  apply  Theorem  3.4 
to  the  sequence  {T„,  n  >  1}  with  the  choice  of  an  =  n  and  mn  ~  0.  Notice  that  7n(0)  = 
— 0(ro),  where  is  defined  by  (4-3).  Thus  from  the  conclusion  (3-23)  of  Theorem  3.4 
and  (4-1)  we  get 


P(*n  >  r0) 


exp(n0(ro)) 


(4-10) 


(2tt  n)fc/2  nLi  to,|V^(t0)|1/2 

for  any  r0  satisfying  (4-5),  where  ^>,  r0  and  |V2^(r0)|  are  given  by  (4-3),  (4-4)  and  (4-6) 
respectively. 


Another  application  of  Theorem  3.4  can  be  given  for  tail  probabilities  of  a  multi¬ 
nomial  distribution.  We  do  not  present  the  details  of  this  application  since  it  is  very 
similar  to  the  above  application  to  the  multivariate  F-distribution. 
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